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THE HERCULES SATELLITE: A STELLAR STREAM IN THE MILKY WAY HALO? 
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ABSTRACT 

We investigate the possibility that the recently discovered Hercules Milky Way satellite is in fact a 
stellar stream in formation, thereby explaining its very elongated shape with an axis ratio of 3 to 1. 
Under the assumption that Hercules is a stellar stream and that its stars are flowing along the orbit of 
its progenitor, we find an orbit that would have recently brought the system close enough to the Milky 
Way to induce its disruption and transformation from a bound dwarf galaxy into a stellar stream. 
The application of simple analytical techniques to the tentative radial velocity gradient observed in 
the satellite provides tight constraints on the tangential velocity of the system (v t = —16^2 kms -1 in 
the Galactic Standard of Rest). Combined with its large receding velocity, the determined tangential 
velocity yields an orbit with a small pericentric distance (i? pG ri = 6^2 kpc). Tidal disruption is 
therefore a valid scenario for explaining the extreme shape of Hercules. The increase in the mean 
flattening of dwarf galaxies as one considers fainter systems could therefore be the impact of a few of 
these satellites not being bound stellar systems dominated by dark matter but, in fact, stellar streams 
in formation, shedding their stars in the Milky Way's stellar halo. 

Subject headings: Local Group — galaxies: dwarf 



1. INTRODUCTION 

The recent discoveries of numerous f aint stellar systems 
around the Milky Way (MW; see e.g. iMartin et atll2008l 
and references therein) has thoroughly changed our view 
of its satellite system. The relatively high velocity disper- 
sions measured from individual stars in these objects are 
usually seen as the sign that they are h i ghly dark mat- 
ter d ominated (e.g. IMartin et al.l I2007L iSimon & Geha 
l2007f l and that they could play a significant role in ex- 
plaining the apparent discrepancy between the number 
of dark matter subhalos seen in ACDM simulations com- 
pared to obser ved luminous dwarf galaxies in our sur- 
roundings fe.g.lTollerud et aJJl2008t IKoposov et~aTll2009t 
iMaccio et alJl201fl ). 

In spite of an ever-increasing knowledge of the prop- 
erties of these newly discovered systems, their origin 
and nature still remain subject to interpretation. Are 
they 'simply' faint equivalents of the previously known 
dwarf galaxies such as Draco, Sculptor or Fornax, lead- 
ing to crucial constraints on galaxy formation and evo- 
lution? Are they disrupted/disrupting versions of for- 
merly brighter progenitors that suffered a destructive 
fate from tidal interactions with the Milky Way? Or do 
they constitute a population of previously unknown sys- 
tems, shaped by different formation and evolution mech- 
anisms, that has remained hidden until the advent of 
systematic searches based on large surveys of the night 
sky? 

The last scenario appears the least likely, given that 
the properties of the faint satellites are a continuous ex- 
tension of our knowledge of Local Group dwarf galax- 
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ies: their sizes and central surface brightnesses are com- 
parable j_although_Jhey_also extend to smaller /fainter 
scales (|Martin et al.l [2008); they appear to globally 
follow the metal licity- luminosity r elation followed by 
brighter galaxies (jKirbv et alJl2008D; they appear to in- 
habit similar dark m atter halos (e.g lStrigari et aT1l2008l : 
I Walker et all 12009). Consequently, it seems more nat- 
ural to envision them as extremely faint dwarf galax- 
ies, as dwarf galaxy remnants, or as a combination of 
both. However, if one is to assume that all recently dis - 
covered systems are dwarf galaxies, IMartin et al.l (2008) 
have shown that the faint MW satellite dwarf galaxies 
are significantly flatter than the brighter ones (observed 
mean ellipticity of (e) = 0.47±0.03 vs. (e) = 0.32±0.02). 

In a thorough analysis of the observational conse- 
quences o f the tidal int e ractio n of a dwarf galaxy with 
the MW, IMufioz et all (|2008f ) show that there is only 
a transient increase in the ellipticity of the satellite as 
it interacts with its host. They nevertheless show that, 
when the satellite is in the final throes of its destruc- 
tion, it can show ellipticities as high as ~ 0.7. In this 
case, the system has lost more than 90% of its stars and 
is becoming unbound, thereby dissolving into a stellar 
stream. In this context, it is interesting to note that the 
large ellipticity measured for systems with My > —7.5 is 
mainly driven by three satellites (about a quarter of the 
sample) that are among the most flattened of all Galac- 
tic satellites: Ursa Major I (UMal; e = 0.80 ± 0.04), 
Hercules (e = 0.68±g;g|) and Ursa Major II (UMall; 
e = 0.63^0 of). These values are reminiscent of those 
measured bv IMufioz et a l. (2008) and suggest that these 
three systems could be transforming into stellar streams 
after their last, destructive, pericentric passage that 
brought them too close to the MW. 

The morphology of these systems is also unlike what 
is seen for rounder dwarf galaxies: deep photomet- 
ric follow-up observations confirm them all to be very 
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elongated and with somewhat distorted morphologies 
(iColeman et alj|2007t lOkamoto et aTll2008l iMufioz et all 
120101 although see iMartin et all 120081 for the large im- 
pact of noise on their distorted shape). This is consistent 
with ste l lar sy stems becoming unbound, as studied by 
IKrounal (fl997l fl. UMall even shows a power-law radial 
dens ity profile, which is typical of a disrupting s ystem 
(e.g. I Johnston et allll~99l iPeharrubia et all [20091 Re- 
gardless of their shape and structure, the possibility that 
these systems are transforming into streams is also tied to 
them having a pericentric distance that is small enough 
for them to be destroyed by the MW's tidal forces. If this 
could easily be the case for UMall that is currently a t 
a heliocentric distance of ~ 30kpc (Zuc ker et all 120061) . 
the viability of an orbit with a small pericenter needs 
to be investigated further for U Mal and Hercules, re- 
siding at distances of 97 ± 4kpc (Okam oto et all [20081) 
and 138 ± 7 kpc (from an error -wei ghted averaging o f the 
values from lAden et al.ll2009al and lSand et al.ll2009D . re- 
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spectively. UMal unfortunately has a small radial ve 
locity with respect to the Galactic Standard of Res 
(y r ~ -lOkms-^IKlevna et al.|[200 l IMartin et~al| [200 
ISimon fc Gehall2007| ). which places it close to its pericen- 
ter or apocenter and limits the constraints that one can 
place on its orbit. Hercules, on the ot her hand, has a 
large receding velocity (y T = 145 kms" 1 ; ISimon fc Gehal 
120071 1 Aden et aLll2009bl ) that could be the consequence 
of a very radial orbit. 

In this paper, we solve for the orbit of Hercules with 
the assumption that it is a stellar stream, in other words 
a disrupting dwarf galaxy that is no longer bound. Our 
goal is to verify that, under these conditions, there is in- 
deed a viable orbit for the system; that is, an orbit that 
brings it close to the MW and can therefore explain its 
transformation into a stream. We use simple dynamical 
arguments based on the observatio n of a tentative ra- 
dial velocity gradient in the system ([Aden et all [2009b) 
to show that there is only a restricted range of tangen- 
tial velocities — and consequently orbits — allowed for 
Hercules if it is a stellar stream. The pericenter of Her- 
cules is very small, thereby confirming that this satellite 
could well be an unbound stellar system. The paper is 
organized as follows: in Section 2, we determine the prop- 
erties of the orbit of Hercules under the hypothesis that 
it is a stream, Section 3 discusses our findings and inves- 
tigates possible discrepancies between a Hercules stellar 
stream and current observations of the system, while Sec- 
tion 4 concludes this work. 

2. DETERMINING THE ORBIT OF THE 
HERCULES STREAM 

The basic assumption of this paper is that Hercules is 
no longer a dwarf galaxy, but instead a disrupted stel- 
lar system that is transforming into a stream. In this 
case, the system has clearly been stripped of all of its 
dark matter and its stars do not form a bound system 
anymore but are instead freely streaming along the or- 

4 See, for instance, the comparison of th e map of Kroup a 's sim - 
ulation with t hat of Hercules obtai ned by Coleman et alj 120071 ). 
as presented inpKroupa et al. ( 20ToT) . 

5 Unless specified otherwise, we use Galactic Standard of Rest 
velocities throughout this paper; that is, velocities observed from 
the position of t he Sun but corrected fo r the motion of the Sun 
around the MW UDehnen fc Binneviri998r) . 
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Fig. 1. — Sketch of the Hercules system in the plane of the orbit 
(left) and in the plane of the sky (right). The variables we use in 
this paper are indicated: the velocity vector v, decomposed into 
the radial velocity, v T , and the tangential velocity of the system, 
v t ; the distance to the center of Hercules, D\ the angular distance 
along the orbit, measured from the center of Hercules, \\ and its 
position angle, 9, defined from equatorial North to East. In the 
plane of the sky, North is to the top and East to the left. 



bit of its progenitoiQ. One thus expects a radial ve- 
locity gradient along the maj or axis of the system that 
Uin fc Lvnden-Belll (j2007l ) andE] (l200cl have form alized 
in the case of generic orbits (see also iBinne vl l2"008l ). We 
refer the reader to these aforementioned sources for more 
detail but, in summary, if we are to follow the naming 
conventions shown in the sketch of Figure Q] and denote \ 
to be the angle along the orbit (measured from the cen- 
ter of Hercules and increasing towards decreasing right 
ascension), D the heliocentric distance to Hercules, v r 
and Vt the radial and tangential velocities of the orbit 
in the Galactic Standard of Rest, and V r -0 the gradient 
of the potential along the line of sight, then the radial 
velocity gradient along Hercules, dw r /dx, is a solution to 
the following quadratic equation: 



^ = Ut + (V r ^. 



(1) 



This results in the following two tangential velocity so- 
lutions of the orbit: 



Vt= 2 



1 / dv T 



(2) 



Since the distance D to Hercules is known and the 
radial gradient of the potential V r "0 can be modeled, a 
direct measurement of a velocity gradient over the body 
of the satellite translates to only two possible orbits, cor- 
responding to the positive and negative solutions of i>t . 

6 Although the stars of a stripped system do not exactly follow 
the orbit of the progenitor, this simplification has little impact on 
our analysis. 
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2.1. Modeling the velocity distribution of Hercules 

stars 

With this direct relationship between the radial veloc- 
ity gradient and the orbit of Hercules in mind, we revisit 
th e determinat i on of a radial velocity gradient presented 
by lAden et all (|2009bl ). From a sample of 18 carefully 
selected Hercules member stars, these authors show the 
tentative presence of a radial velocity gradient in a direc- 
tion roughly consistent with the major axis of the system. 

2.1.1. The model 

We employ a maximum likelihood algorithm to fit a ve- 
locity gradient through the same data points they have 
used. Our model has four parameters: the radial ve- 
locity gradient, du r /dx, the mean radial velocity at the 
center of Hercules tv, the direction of the linear gradient 
on the sky, defined by its position angle from equatorial 
North to East, 9, and the velocity dispersion, s, of mem- 
ber stars around the velocity gradient. The goal of the 
maximum likelihood technique is to find the set of these 
four parameters that maximizes the likelihood function 

( 3 ) 

where li (di> r /dx,I^, 9, s) is the probability of finding the 
datum i given the set of parameters. In the current 
problem, each star is defined by its radial velocity, v tt i, 
and associated uncertainty, v err ^, as well as its right as- 
cension and declination (ai,Si) that we convert to its 
distance from the center of Hercules (ao,<5o) such that 
Xi = {cti — ao)cos((5o) and Yi — <5j — Sq. The angular 
distance of star i along an axis of position angle 9 is then 
Di = XiSm(9) + YiCOs(9), yielding a difference between 
the modeled and measured radial velocities at this posi- 
tion of 

Au r ,i = v I:i - (j^Vi + ^) ■ ( 4 ) 

This velocity difference around the radial velocity gradi- 
ent is finally modeled by a Gaussian, whose standard de- 
viation is the velocity dispersion of the system, s, added 
in quadrature to the uncertainty in the velocity measure- 
ment of star i. This leads to the following expression for 
if 

li {iT^^'*) = / 1 = CX P \-\ 2 i f 2 I ' 

Now that the likelihood function is entirely expressed 
as functions of the parameters and the data point prop- 
erties, the model that maximizes it is determined by ex- 
ploring a fine grid over the four dimensions of the param- 
eter spac^H- In order to determine the uncertainties on 
the measurement of a given parameter pj , we marginal- 
ize over the other three parameters and assume that the 

7 One could wonder whether a model with a radial velocity gra- 
dient that has one more parameters than a model with a fixed 
mean velocity and velocity dispersion is warranted when only 18 
stars are available. A likelihood-ratio test in fact yields that the 
simpler model can be rejected at the 94% confidence-level. 
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Fig. 2. — Two-dimensional likelihood contours for parameters of 
our model. The filled circle represents the best model for each 
parameter set. From this point outwards, the contours indicate 
drops in likelihood of 50%, 90% and 99%. These are very regular, 
showing that the maxima of the likelihood distributions are well 
defined. 



marginalized likelihood function C is well-behaved near 
the best model (i.e, Gaussian-like or close to Gaussian). 
This allows us to use the property that 21n(£') behaves 
as a x 2 distribution and determine the k-a confidence 
interval of parameter pj as being bound by the values 
of pj that correspond to 21n(£') dropping by k 2 . The 
uncertainties given below correspond to this definition of 
the Ter confidence interval. 

2.1.2. Results 

We have checked that, if we are to force no velocity dis- 
persion for t he mo del (s = Okms -1 ), the best model of 
lAden et al.1 ((2009b) is recovered with a heliocentric «7 = 
45.li0.4kms" 1 , and dv r /dx = 15.4 ± 2.8 kms" 1 kpc" 1 
for a position angle of 9 — —35°. If all the parameters 
are allowed to evolve freely, they converge on the fol- 
lowing best values: heliocentric vT- = 45.0 ± 1.1 km s -1 , 
s = 3.5to;gkms _1 and dv r /dx = 14.1^7 jj km s _1 kpc -1 
(as part of the velocity gradient is accounted for by the 
assumed internal velocity dispersion), for a position an- 
gle of 9 = -37°. 

However, under our assumption that Hercules is a stel- 
lar stream, the axis of the radial velocity gradient must 
be the axis along which the system is elongated, that 
is, its major axis. Conseq uently, the position angle of 
the system, 9 = -78 ± 4° ([Martin et al.ll2008H . is used 
as a prior in order to find the best model, whose two- 
dimensional marginalized likelihood functions are shown 
in Figure [2] for the three remaining parameters. From 
the one-dimensional marginalized likelihood functions, 
we determine the following values for the best model: 

= 144.7il.2kms" 1 , s = 3.5±J;g kms" 1 and dv r /dx = 
10.2 ±6.0 km s" 1 kpc- 1 = 24.6 ± 14.5 kms" 1 deg" 1 . 

2.2. The distribution of likely Hercules orbits 
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dv r /dx (km/s/deg) v t (km/s) i? peri (kpc) 

Fig. 3. — Marginalized relative likelihood distributions for the radial velocity gradient, tangential velocity and perigalacticon, whose best 
values are dti r /dx = 25 ± 15kms — 1 deg -1 , vt = — 16^2 kms - 1 and -R per i = kpc, respectively. From top to bottom, the dashed lines 
intersect the likelihood functions at the boundaries of the 1-, 2- and 3-cr confidence intervals. 



Having obtained the likelihood distribution of the ve- 
locity gradient, shown in its marginalized form in the left 
panel of Figure [3J it is now possible to use equation ([2]) 
to determine the likelihood distribution of the tangential 
velocity of the Hercules orbit. The distance is assumed to 
be D — 138 ± 7 k pc, obtained b y avera ging t he distance 
measu rements of lAden et all (|2009al ) and iSand et all 
(2009). The Milky Way potential we place ourselves in 
is a mixture model whic h combines the M iyamoto-Nagai 
disk and bulge defined bv lPaczvhskil (|1990f) wit h the adia- 
batically contracted NFW halo constrained bv iXue et al.l 
(|2008h . leadingto V r <0 = -170(kms- 1 ) 2 kpc -1 at the lo- 
cation of Hercules. Choosing the positive solution for i\ 
yields very large velocities that produce unbound orbits 
for Hercules with a very large pericenter. This is clearly 
not a viable solution for a system disrupted by a recent 
passage close to the Milky Way center. The negative so- 
lution for vt , on the other hand, yields much more inter- 
esting results. The marginalized likelihood distribution 
for this case is shown in the middle panel of Figure [3] and 
reveals a preferred velocity of Vt = — I6I22 kms -1 , con- 
sistent with a very radial orbit given the mean velocity 
of Hercules (v^ — 144.7 ± 1.2kms _1 ) at its center. 

Given that the sky position, radial velocity, distance to 
and direction of motion of the assumed Hercules stream 
are known, determining the tangential velocity likelihood 
distribution of the system is equivalent to knowing its or- 
bit likelihood distribution. In particular, there is a direct 
relationship between C(v t ) and the likelihood distribu- 
tion of the stream's Galactocentric distance at the last 
pericenter, £'(i? por i). By integrating the orbits for 2Gyr 
backward and forward in time in the potential described 
above, we obtain the distribution of £'(i? pcr i), also shown 
in Figure 02 from which we determine tight constraints 
on the small pericenter of the orbit: i? pol -i = 6^2 kpc. 
This pericenter was reached ~ 0.6 Gyr ago. The prop- 
erties of the corresponding preferred orbit are shown in 
Figure HJ 

3. DISCUSSION 

Under the assumption that the radial velocity gradient 
observed along the body of Hercules is due to the sys- 
tem being a stellar stream, we are able to find an orbit 
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Fig. 4. — Orbit of the Hercules stellar stream using our best 
model as given in the text, integrated forward and backward from 
the present location for 2 Gyr in each direction. From top to bot- 
tom, the panels show the location of the orbit on the sky in Galactic 
coordinates, as well as its heliocentric distance and radial velocity 
relative to the Galactic Standard of Rest. The solid and dotted 
lines respectively denote the forward and backward-integrated or- 
bit from the present location of Hercules, which is indicated by an 
arrowhead in the direction of the orbit. 



that has brought it to within 6_2 kpc of the Galactic cen- 
ter. But is this close enough to induce the destruction 
of a dwarf galaxy into a stellar stream during the last 
0.6 Gyr and does the very radial orbit, with a perigalac- 
ticon to apogalacticon ratio of 6:219 (i.e. an eccentricity 
e = 0.95), contradict the current knowledge of Hercules? 
We now discuss these two points and the various conse- 
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quences of the orbit we have derived above on the obser- 
vational properties of the stellar system. 

3.1. A small enough pericenter for tidal disruption? 

Whether or not a pericenter of 6^2 kpc is small enough 
to induce the tidal destruction of Hercules is a function 
of the properties of the stellar system at that time of its 
journey in the Milky Way potential and, as such, out of 
reach to the current observer. We can however use ob- 
served properties of Hercules to conservatively estimate 
the minimum pericentric radius it should have reached to 
undergo a strong tidal intera ction with the Milky Way. 



By definition (jKingl ITT)li2i. at its pericenter R 1 
satellite is pruned to its tidal radius, rt, such that 



pern 



n ~ R 



peri 



M(i? pcri )(3 + e) 



1/3 



(G) 



where e is the eccentricity of the satellite's orbit around 
the Galaxy, m is the total mass of the satellite and 
M (i?p 0r i) is the mass of the Milky Way within the peri- 
centric radius. Assuming a constant circular velocity v c 
(expressed in kms^ 1 ) for the plausible range of Hercules 
pericentric distances, this mass can be expressed (with 
i?peri expressed in kpc) as 



Af(i? pori )^%^xl0 6 M G 



(7) 



Final ly, if we assume that v c ~ 200 km s -1 (X ue et all 
120081) . inserting equation (JT]) into equation © and solve 
for -Rpcri, we obtain 



i? pcri - r\l 2 ( £±£> x 10- 



1/2 



(8) 



with m expressed in solar units. 

We are left with a function of the properties of Her- 
cules at pericenter and we know that, for the favored 
orbit determined above, e = 0.95. The two remaining 
unknowns need to be assumed. In the case of the tidal 
radius at pericenter, advocating for tidal disruption re- 
quires it to be much smaller than the current extent of 
the syst em that has an o n-sky King limiting radius of 
1.4 kpc (|Sand et al.l 12009ft . We assume r t ~ 0.3kpc, a 
value that is close to the current half-light r adius of Her- 
cules (jMartin et all 120081: ISand et al.ll2009D . as a neces- 
sary condition for the system to be stripped of a large 
fraction of its stars and transformed into stellar stream. 
Regarding t he ma ss, we use the value determined by 
lAden et al.l (|2009b) within 300 pc: m ~ 2 x 10 6 M Q . This 
leads to R 



pen 



23 kpc being required for a strong inter- 
action to take place between Hercules and the Milky Way. 
This value is of course very uncertain as it is strongly de- 
pendent on our choice of r t and m but it confirms that, 
with i?p C ri = 6^2 kpc, the preferred orbit determined in 



Section 12.21 appears to have a small enough pericenter 
to induce the transformation of Hercules into a stellar 
stream. 

3.2. The shape of Hercules 

With the distance to Hercules being much larger than 
the Sun-to-Galactic-center distance, the system's or- 



bit being very radial means that its physical, three- 
dimensional size is much larger than its projected size 
as measured on the plane of sky. Transforming the mea- 
sured half-light radius and ellipticity measured on the 
sky, r/i and e = 1 — b/a, into their de-projected equiva- 
lents, r^deproj and edeproj, can be done easily from simple 
geometrical considerations obtained directly from Fig- 
ure [1] 



JXdeproj = Hi (COS \) = r h 1 + -§■ 



1/2 



^deproj 1 I 



-1/2 



(9) 
(10) 



From on-sky measurements ( r^ = 230 ± 30 pc and 
e = 0.67 ± 0.03; ISand et alJfeOOa ). these relations yield a 
de-projected half-light radius and ellipticity of r^deproj = 
1.5 kpc and edeproj = 0.95. One should note, however, 
that in the case of a stellar stream, the half-light radius 
is not a very meaningful quantity as the stream's surface 
brightness profile is not expected to follow the typical 
exponential, Plummer or King profiles that have so far 
been used to derive the properties of this satellite. Dis- 
rupting systems are expected to f ollow power-law den- 
sity profiles in their ou tskirts (e.g. Uohnston et al.lll999l : 
iPenarrubia et al J 12009). but current data have not yet 
been used to quantify the validity of such a density model 
for Hercules. 

3.3. The issue of the distance gradient 

If Hercules is a stellar stream as opposed to a dwarf 
galaxy, its stars are, as mentioned before, freely stream- 
ing along its orbit. This has the consequence of naturally 
producing the radial velocity gradient that was used in 
Section 2 to constrain the orbit, but it also produces a 
distance gradient along the body of the system that is 
all the more i mport ant as the orbit is very radial. In 
Uin fc Martini (|200ll . we have studied this distance gra- 
dient in the generic case and shown that, following the 
notations used here, it can be expressed as: 



dx v t ' 



(11) 



With the properties of the Hercules orbit determined 
above, it follows that the system should exhibit a dis- 
tance gradient of dD/dx = — 22l3 kpcdeg~ 1 . 

The current best estimate of a distance gradient, based 
on color-magnitude diagram (CMD) fitting of small re- 
gions on e ither side of Hercu les' center, has been deter- 
mined by ISand et al.l (|2009f >- They find that the dis- 
tance gradient over the system is at most of ~ 6 kpc over 
~ 6arcmin, or dD/dx — GOkpcdeg" 1 . This gradient 
goes in the opposite direction t o the o ne we determine 
from our orbit, but ISand et "ail (|2009f ) show that their 
detection of a distance gradient is very tentative. In 
fact, from their Figures 13 and 14, one can see that it 
is just as likely for an inverse gradient to be present in 
their data. Thus, there is no outstanding discrepancy 
between the current observations of the system and the 
distance gradient that should be observed if Hercules fol- 
lows the orbit we propose in this paper. However, given 
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the low density of Hercules' stars, it will be difficult to 
observe the predicted dis tance gradient. Co ntrary to, 
for instance, Draco (e.g. iKlessen et al.l l2003h . the sys- 
tem hosts but a few horizontal branch stars from which 
an accurate distance gradient measurement could be de- 
rived and the alternative of using deep data reaching the 
main sequence turn-off is plague d by contaminati on at 
the faint end of the CMD (cf. the ISand et~al"1l2009l anal- 
ysis). Even though the orbit of the Hercules stream is 
very radial, confirming the presence of its large distance 
gradient will be a trying task. A search for the few RR 
Lyrae stars that should inhabit the system might be the 
best alley of investigation to ever constrain the presence 
of a distance gradient. 

3.4. A tidally disrupted Hercules? 

If there is no evidence that Hercules cannot be a stream 
in formation, there are in fact some signs that it could not 
be a pristine, bound s tellar system. The deepest avail- 
able data of Hercules (jColeman et al.l l2007t iSand et all 
2009) show the presence of some material along the ma- 
jor axi s of the system, e xtending away from its main 
body. ISand et"al1 (|2009t l highlight that this extension 
has a CMD that is tentatively similar to that of Hercules 
and also coincides with a small overdensity of likely Her- 
cules blue horizontal-branch stars. In addition, further 
along the major axis (or likely orbit of the stream), they 
investigate the nature of three small clumps of stars that 
appear in their smoothed maps within ~ 0.5° of Her- 
cules' centroid, and again find tentative evidence of a 
c onnection to the ste llar system. 

lAden et al.l (|2009bD also find some asymmetry in the 
distribution of their carefully selected Hercules member 
red giant branch stars that they ascribe to the possible 
effect of tides. Although this extension of three stars 
is not aligned with the major axis and therefore seems 
discrepant with our scenario of Hercules being a stellar 
stream, one has to note that their photometric analysis 
only covers a limited field of view that does not favor the 
detection of stars extending along the major axis of the 
system. Thus, it is not inconsistent with the scenario we 
propose. 

In addition, and as mentioned in the introduction, Her- 
cules has a st riking resemblance with the result of the 
Kroupa (1997) simulations which emphasize that, even 
though they are not systems in equilibrium at the center 
of massive dark matter halos, disrupting stellar systems 
can have properties consistent with those observed for 
dwarf galaxies. A rece nt, more detaile d re-si mulation 
of this initial work by iMetz fc Kroupal (|2007l ) further 
shows that such disrupting systems can share the ob- 
served properties of the population of recently discov- 
ered, faint, MW dwarf galaxies. There are therefore ways 
to form a disrupting stellar system that shares the ob- 
served properties of Hercules. 

This interpretation seems to contradict that of 
iPenarrubia et al.l ([2008), who find that the recently dis- 
covered satellites, including Hercules, are unlikely to be 
the remnants of brighter MW dwarf galaxies that once 
shared the current properties of Draco, Fornax, or Sagit- 
tarius. Their conclusion mainly stems from the relatively 
large velocit y dispersions of the newly discovered stel- 
lar systems. lAden et al.l ()2009al ) have, however, shown 
that Hercules in fact has a low velocity dispersion that 



we find to be even smaller when one properly accounts 
for the observed velocity gradient in the system. Fol- 
lowing b ackwards the 'tidal evo lutionary tracks' mea- 
sured by fPenarrub ia et al.l (|2008l ) to track the evolution 
of the properties of observed dwarf galaxies embedded in 
ACDM dark matter halos as they undergo tidal disrup- 
tion (their Figure 10), the current velocity dispersion and 
luminosity of Hercules (s,L) = (3.5 kms -1 , 3.6 x 10 4 L Q ) 
correspond to the properties expected for an initial satel- 
lite with (s,L) = (8.7kms _1 ,3.6x I0 5 L Q ) that has been 
stripped of 90% of its stellar masfl These properties are 
in fact very similar to the current velocity dispersion and 
luminosity of the Draco or Ursa Minor dwarf galaxies, 
making it possible that Hercules is actually the remnant 
of such a galaxy. 

However, if this is truly what is happening to Her- 
cules, the low metallicity measured for its stars (mean 
[Fe/H] between —2.6 and —2.0, depending on the study ; 
Simon fc Ge ha 2007; Kir bv et al.l feOOS: Ko ch et al.l feOOS: 
Aden et al.ll2009aft and its broad agreement with the lu- 
minosity / metallici ty relation followed by dwarf galaxies 
(jKirbv et al.l[200l would seem to imply that the progen- 
itor of the observed system could not have been much 
brigh ter than the current luminosity of He rcules (My ~ 
-6.5; lMartin et al.ll2008HSand et al.l(2009l ). In fact, the 
large spread in the relation would nevertheless make it 
possible for the observed stellar system to be a disrupted 
versio n of Draco or Ursa Minor ([Fe/H] ~ — 2.0: IWinnickl 
120031 ). A violent tidal interaction with the Milky Way, 
which is consistent with the very eccentric orbit we have 
determined (e = 0.95), would also allow for a system 
which quickly goes from bound to disrupted, without the 
need for a long, gradual peeling off of its stars that would 
allow the progenitor to stay bound for a longer period of 
time. 

There remains the issue of timing and how likely it is 
to be observing a system that is exactly in its disrup- 
tion phase but that is still concentrated enough for it to 
be observed as an overdensity of stars. This could make 
for an uncomfortabl e coincidence but, a s has been previ- 
ously mentioned bv lAden et al.l (j2009b[ ). rFcllhaucr et al.l 
(2007) have shown in simulations of the UMall stellar 
system that this process requires ~ 1 Gyr. With the 
last, destructive pericenter of the orbit we have deter- 
mined happening only ~ 0.6 Gyr ago, it would therefore 
not be surprising that we could still be observing the 
core of a Hercules system currently transforming into a 
stream. The UMall simulation with its different orbit is 
of course not directly comparable to the Hercules case, 
and detailed simulations of the evolution of a disrupting 
Hercules following the orbit we have determined are war- 
ranted, but we currently find no outstanding issues with 
the scenario taken as the assumption of this paper. 

3.5. Hercules mass estimates 

By removing foreground Milky Way stars that happen 
to share velocities similar to that of member stars, the 
velocity dispersion of Hercules and hence its mass have 
been severely corrected downwards to yield a dynami- 
cal mass within 300 pc of only M 300 = 1.9±oll x 10 6 M 

8 This ~ 90% loss in stellar mass is a nece s sary c ondition for 
the dwarf galaxies simulated by Muhoz et al. (2008) to become 
unbound and show large ellipticities. 
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(|Aden et al.l 2009bf). This value is low enough that it 



falls off the iStrigari et al.l (|2008l ) 'common mass-scale' 
of dwarf galaxies that all share M300 ~ 1O 7 M0. This 
scale can easily be explained by dwar f galaxy forma- 
tion models (e.g. iMaccio et all 120081: iLi et all 120091: 
lOkamoto fc Frenkll2009D and the Hercules dwarf galaxy 
being a significant outlier could be worrying. 

One must nevertheless remember that there are nu- 
merous assumptions that enter the determination of such 
dynamical masses, first and foremost that the stellar sys- 
tem is in virial equilibrium. This is obviously not a valid 
assumption if the system is currently disrupting and this 
would lead to spurious mass estimates. Taken in conjunc- 
tion with the tentative observation of a velocity gradient 
and the presence of some stars possibly streaming out 
of the system, the abnormally low velocity dispersion 
of Hercules could be yet another sign that it is in fact 
a stream in formation, as opposed to a peculiar bound 
dwarf galaxy. 

4. CONCLUSION 

We have shown that, under the assumption that the 
tentative velocity gradient observed in Hercules is a sign 
that it is no longer a bound dwarf galaxy but instead 
a stellar stream in formation, there is a viable orbit 
that can explain the properties of the system. We con- 
strain the tangential velocity of Hercules to be only 
v t = —16^22 km s _1 , which makes for a very eccen- 
tric orbit (e = 0.95) whose pericentric distance is only 
-Rpcri = 6^2 kpc. This value is small enough that it could 
explain the tidal disruption of a satellite that we cur- 
rently see as a stellar stream in formation after its last 
destructive pericentric passage close to the MW center. 

We do not claim that there can be no other scenario 
that could explain the observed properties of the system 



and the hypothesis of a stream in formation that we have 
presented in this paper may not necessarily be the plight 
of Hercules. Our intention is to provide avenues of inves- 
tigation in order to understand the peculiar properties of 
some of the recently discovered faint dwarf galaxies. We 
show that, in the case of Hercules, tidal destruction seems 
to be a perfectly viable option to explain its surprisingly 
large ellipticity and low velocity dispersion. It is worth 
emphasizing that we are not advocating that all faint 
dwarf galaxies are stellar streams but simply that some 
systems, namely those with very large ellipticities (e.g. 
Hercules, UMal, UMall), could well be shaped by such 
a process. This hypothesis can, and should, be tested by 
deeper and wider photometric follow-up that can reveal 
thes e systems to be muc h larger than originally measured 
(e.g. IMunoz et al.1 12010 for UMall) and/or, in the case of 
Hercules, by more extended spectroscopic coverage that 
would definitely confirm (or rule out) the presence of the 
radial velocity gradient that is our core assumption. 

If the systems with the largest measured ellipticities 
are indeed stellar streams in formation, that would alle- 
viate the necessity of having to find a galaxy formation 
mechanism that produces more flattened systems at the 
faint end of the luminosity function. This simple possi- 
bility that Hercules could be an unbound system should 
also make one wary of treating the population of newly 
discovered stellar systems as a single population of ob- 
jects. 
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